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I. INTRODUCTION 



The lineshape scan of the threshold top pair production cross section constitutes an 
integral part of the top quark physics program at a future e~^e~ or 77 collider. Q, |^ 
Because in the Standard Model the top quark width Ff ^ 1.5 GeV is larger than the typical 
hadronization energy Aqcd? it is expected that the lineshape is a smooth function of the 
c. m. energy, and that non-perturbative effects are strongly suppressed. From the rise of 
the cross section a precise measurement of the top quark mass will be possible, while from 
the shape and the normalization of the cross section one can extract the top quark Yukawa 
coupling yt, the top width and the strong coupling. \^ Past fixed order next-to-next-to- 
leading order (NNLO) computations of the cross section have shown that a measurement of 
the top quark mass in a threshold mass scheme with theoretical uncertainties of 100-200 MeV 
or better are feasible. However, in the fixed order approach the theoretical uncertainty 
of the normalization of the NNLO cross section were estimated at the 20% level, which 
would jeopardize competitive measurements of yt, F^ or 

Recently, the renormalization-group-improved e~^e~ top threshold cross section was com- 
puted i^ the framework of an effective theory for non-relativistic heavy quark pairs, 
called vNRQCD. & i H H E Q This effective theory describes the dynamics of heavy 



quarkonium systems, when the hierarchy of scales m ^ mv ^ mv"^ ^ Aqcd is satisfied, m 
being the mass and v the average cm. velocity of the quarks. The matching is carried out 
at the hard scale fi = m onto a potential-like theory with both soft and ultrasoft degrees of 
freedom. For scales fi < m the correlation of energy and momenta is accounted for since the 
ultrasoft and soft renormalization scales fiu and fis are related, fiu = fi'^/m = mu'^. The 
running in vNRQCD is expressed in the dimensionless scaling parameter u. All operators 
(and their coefficients) are evolved from z/ = 1 to z/ ~ v of order of the average cm. velocity 
of the quarks, where matrix elements are free from large logarithmic terms. In dimensional 
regularization the factors of fil; and /i^ multiplying operators in the renormalized effective 
Lagrangian are determined uniquely from the v power counting in d dimensions. 0| 

In renormalization-group-improved perturbation theory the expansion of the normalized 
cross section R takes the parametric form 

R = = (— )^^(a,lnw)' X |l (LL);a„w(NLL);a^^a,f,f2(NNLL)| , (1) 

where f <^ 1 is the top quark velocity and where the indicated terms are of leading log- 
arithmic (LL), next-to-leading logarithmic (NLL), and next-to- next-to-leading logarithmic 
(NNLL) order. In Refs. U all logarithms were summed in the Wilson coefficients 
of the operators that contribute to the cross section at NNLL order except for the Wilson 
coefficient ci of the leading order spin-triplet current 

Ji,p = '^(^^2)x-p , (2) 

for which only the NLL anomalous dimension was known. §]. In Refs. jol, 0, 0| it was 
shown that the summation of logarithms leads to a significant reduction of the normalization 
uncertainties, and a theoretical uncertainty of 3% was estimated for the NNLL cross section 
in e+e~ annihilation. The computation of the NNLL anomalous dimension of the current 
Ji^p, is an important task for the determination of the full NNLL order cross section and 
for a cross check of the error estimate made in Refs. [gI. ItI. Es^ . 
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The evolution of ci(z/) is obtained by integrating the anomalous dimension 



d 



NNLL^ 

ci 



The LL order anomalous dimension is zero. The NLL order term reads 



(3) 



7, 



,NLL/ 
ci \ 



167r2 



(4) 



where = 2 is the squared quark total spin operator for the spin-triplet configuration.^ 
The terms Vc**'* and Vg*'' 
of order asV~^ 



V, 



r , V]*'' are the color singlet Wilson coefficients of the 
and ttsV (1/m^, (p^ + p'^)/(2m^k^), S^/m^), respectively 

(1) /o(^) 



le pot 



otentials 
while 



are the Color singlet Wilson coefficients of the sum operators Ol^\ introduced in 



Ref. jlSf]. In the convention for the vNRQCD operator basis used in Refs. j8|,|l3| the potentials 
at order a'^v^ contained operators having the momentum dependence l/|k|, k = p — p' 
being the momentum transfer, which mixed into 7^ at NLL order. For technical reasons 
explained below in Sec. lIIIBl we adopt in this work the convention where all a^v^ potentials 

Instead of the l/|k| potentials 



(1) (T) 

are represented by sum operators in analogy to , 



used in Refs. I8l ll 



and O 



(T) 



AC, 

, respectively. 

Ref. 



giving the contribution 

(1) 



we thus have the sum operators 

l to the vNRQCD Lagrangian with Wilson coefficients V^'^'' and 
At LL order the coefficients are in agreement with the convention in 
13i | and contribute to the second term on the RHS of Eq. (@]). The explicit form for 
the operators is given in the appendix. Note that in the following we frequently refer to 
the sums over field indices in the sum operators as loop integrals in order to simplify the 
presentation. 

The mixing displayed in Eq. (jU arises from two-loop vertex diagrams containing only 
potential loops and insertions of the Coulomb potential, the subleading heavy quark kinetic 



energy operator and the 1/m-suppressed potentials. j8| At NLL order, in Eq. (jU, ag, V, 



'^kik2 "^'fc^'^'' iieed to be known at LL order for all values u < 1. 

HQ 



At NNLL order there are two classes of contributions. The first is due to the two-loop 



mixing shown in Eq. (@]) and requires the NLL results for a^, V, 



•^2 ) 



kl,k2 



and , but it does not modify the form of the NLL anomalous dimension. The second 
class requires the computation of three-loop vertex diagrams with potential loops and either 
soft or ultrasoft loops that require new Ci counterterms. We call this class "non-mixing 
contributions" since it leads to genuinely new contributions in the anomalous dimension of 
Ci. By power counting there are no contributions from diagrams with three potential loops 
or which have both soft and ultrasoft loops. 

In this paper we present the non-mixing contributions to the NNLL anomalous dimen- 
sion of Ci. As mentioned above, all order a^v^ potentials are presented by sum operators. 



In this work we keep the exphcit dependence on the total quark spin operator. Thus the resuhs can be 
generahzed to the spin-singlet configuration where — 0. 
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Otherwise, we use the vNRQCD velocity renormahzation group with the notations and 
conventions of Refs. j^, 0, We employ the MS scheme in d = 4 — 2e dimensions. 

The outline of the paper is as follows: In Sec.|n] we describe the method we have used 
to carry out the computations and illustrate it by rederiving the NLL anomalous dimension 
of C\. Based on renormahzation group invariance we also predict the NNLL order 
term of the renormahzation constant of C\. In Sec. lIIII present and discuss our result for 
the NNLL non-mixing contributions of the anomalous dimension of C\. Section llll Al and 
nil Bl are devoted to the contributions involving the dynamics of ultrasoft and soft degrees 
of freedom, respectively. The contributions to the NNLL renormahzation group equation 
of Ci are determined and discussed in Sec. lIII CI and the modifications of the formulas for 
the NNLL heavy quark pair production cross section at threshold in e^e~ annihilation are 
discussed in Sec. lIII Dl Our conclusions are given in Sec. lIVI The paper has three appendices 
where we have collected formulas for the reader interested in the details of the computations. 



II. METHOD OF THE COMPUTATION 



The standard method to determine the renormahzation constant of the current Ji p con- 
sists of computing the overall UV-divergences of quark-antiquark-to- vacuum on-shell matrix 
elements of spin-triplet currents at a certain loop order including lower order counterterm di- 
agrams needed to subtract the subdivergences. There are, however, technical complications 
in this approach due to the existence of IR-divergent Coulomb phases for on-shell quarks 
and from the fact that the vertex diagrams in general depend on three physical scales, the 
mass m, the c. m. energy E and the external quark momentum p. Note that imposing the 
on-shell condition = mE on the quark-antiquark-to-vacuum matrix elements in dimen- 
sional regularization, yields IR-divergent l/e" poles that are very difficult to separate from 
the UV-divergences. In practice it is therefore necessary to compute the quark-antiquark- 
to-vacuum matrix elements in an asymptotic expansion for p^ — mE <C mE, where the 
IR-divergent Coulomb phases manifest themselves as powers of ln((p^ — mE)/mE). An ef- 
ficient way to avoid these complications and, in addition, to reduce the number of diagrams 
that have to be computed is to consider current correlator graphs rather than the vertex 
diagrams. The correlator graphs are obtained from closing the external quark lines of the 
vertex diagrams with an additional insertion of the current Ji,p. This means that one has to 
determine diagrams with one more loop, but it reduces the number of diagrams to be con- 
sidered, eliminates the quark momentum p as an external scale and avoids the IR-divergent 
Coulomb phases. Moreover, since all IR-divergences cancel, using correlator graphs it is not 
necessary to distinguish on- and off-shell contributions since any off-shell term that could be 
relevant for the renormahzation constant leads to scaleless integrals in dimensional regular- 
ization which automatically vanish. In this approach the three-loop (NNLL) renormahzation 
constant of the current Ji,p is obtained from the subdivergences in four-loop correlator di- 
agrams that remain after the one- and two-loop subdivergences have been subtracted. The 
surviving four-loop overall divergences are canceled by external vacuum-type diagrams and 
not related to the renormahzation of the current Ji,p. The method applies analously at 
any order. Interestingly, the overall divergences of the four-loop correlator graph vanish in 
dimensional regularization because the graphs are non-analytic in the c. m. energy, and there 
are no operators in the effective theory that could absorb the overall divergences. The result 
for the renormahzation constant obtained by this method will agree with the one obtained 
from quark-antiquark-to-vacuum on-shell matrix elements. 



4 



(a) (b) 

ocx> 



(c) 




(d) 



(e) 




FIG. 1: Three-loop current correlator diagrams and the counterterm diagram for the computation 
of the NLL anomalous dimension of ci. 



For illustration let us reconsider the computation of the NLL renormalization constant 
that leads to Eq. The relevant correlator diagrams are shown in Fig.^ Four-quark 
interactions without label refer to the Coulomb potential. Crosses on quark lines refer to 
insertions of the quark kinetic energy at subleading order. Here and throughout the paper, 
combinatorial factors and diagrams obtained by flipping the graphs left-to-right and up-to- 
down are to be understood. There are no one-loop sub divergences that have to be subtracted 
in this case. Adopting the form 



ci + 5ci 



Zci Ci 



(5) 



for the unrenormalized Wilson coefficient of the current Ji^p with the renormalization con- 
stant written as 



NNLL,1 



+ ... , 



(6) 



one obtains 



(7) 



where J^^^ii^) is given in Eq. (@]). Taking into account the velocity renormalization 
group equations for the vNRQCD coefficients in d = 4 — 2e dimensions {v-^gsimi') = 
—egsirni') + . . ., v-^gsimv'^) = — 2e(7s(mz/^) + . . ., etc.) one can derive Eq. (jH) using that 
Co is renormalization group invariant. From the same relation and the fact that all l/e"' 
(n = 1, 2, . . .) terms cancel in the renormalization group equations one also finds the coeffi- 
cient of the term at NNLL order (a^ = Q;s(mz/), a„ = Q;s(mz/^)), 



Sz. 



,NNLL,2 
ci 



24:71 



1927r2 



3-CD + 2ci,(l 



k2 I 



-Cp{ - CF{C2b + 2C; 



2c 



2887r 



Cf < 8C1 + UC'f + CaCf 



f + yc^ + 4(-13 + 7S^) 



(8) 
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FIG. 2: Four-loop graphs for the calculation of the ultrasoft non-mixing contributions of the NNLL 
anomalous dimension of ci. 



where = 11/3C^ — A/STrii is the one-loop QCD beta-function, and and cp are the 
coefficients of the Darwin and the magnetic operators of the HQET action j^. The terms 



^2a,2b,2c 



are coefficients associated to the soft 6-field operators 



(2) 



which were intro- 



duced in Ref. [13[, and which mix into the order a^f potential operators. The formulas for 
the coefficients, which all depend on the scaling parameter u, are collected in App.O Note 
that at NNLL order only the diagrams that lead to non-mixing contributions of the anoma- 
lous dimension of Ci can also contribute to SZ^_^^^''^. Thus, Eq. (jSj) provides a non-trivial 
cross check for our result and the consistency of the effective theory under renormalization. 



III. NON-MIXING CONTRIBUTIONS AT NNLL ORDER 



For the non-mixing contributions to the NNLL anomalous dimension of Ci one can distin- 
guish between four-loop correlator diagrams containing ultrasoft loops and those with soft 
loops. Diagrams with both ultrasoft and soft loops do not exist at this order. Note that 
the ultrasoft and soft non-mixing contributions are separately gauge-invariant since they 
are induced by ultrasoft and soft gluons, respectively, which represent different degrees of 
freedom in the effective theory. 



A. Ultrasoft non-mixing contributions at NNLL order 

For the determination of the ultrasoft non-mixing contributions we adopted the Coulomb 
gauge, where the propagation of the longitudinal ultrasoft gluon field component does 
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(a) 8V2,r,s (b) SV, V2,„ (C) 5Vk 



(d) 6Vt (e) 5V, (f) 8V, 
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FIG. 3: Counterterm graphs for the removal of subdivergences in graphs of Figs. [21 and 0] Graphs 
with wave function renormaUzation constants are to be understood. 



not contribute. For the transverse ultrasoft gluon field A the leading order couplings to 
quarks and the Coulomb potential have to be taken into account. The complete set of four- 
loop correlator diagrams is displayed in Figs. El The diagrams needed to cancel the one-loop 
subdivergences are displayed in Figs.|Hk,c,h using the graphical notation from Refs. P. lol . 
Here, 5c2 denotes the counterterm of the Wilson coefficient of the f ^-suppressed spin-triplet 
current l/rn? t/^^ p^cr(icr2)xlp- SB] There are no two-loop subdivergences. After the removal 
of the one-loop subdivergences the contribution to the renormalization factor from each 
of the diagrams in Figs. El is of the form 



47r 



A B 

- + - 

e 



(9) 



The contributions to A and B as well as the totals are given in Tab.lH The results com- 
ing from graphs |21 and ^Bpv are zero, but non-trivial and have been obtained from explicit 
computation. The results do not involve factors of p = 7^; — ln(47r) because we have imple- 
mented the MS scheme by scaling each term ns or nu with a factor e''^'^. Note that each 
counterterm diagram can in general contribute to several four-loop diagrams. For exam- 
ple, the counterterm diagrams involving 6V2,r,s are required for the subtraction of one-loop 
subdivergences in the four- loop graphs |2KH21i and|2hH21, whereas the diagrams with SVki,k2 
have to be considered for the graphs |2f, I2k and |2k. On the other hand, the diagram with 
Sc2 is needed for the graphs and Et. The graphs Efe? El and Ehi do not have any one- 
or two-loop subdivergences. From the totals given in Tab. HI one finds agreement with the 
prediction of the ultrasoft term made in Eq. (|H1). 

It is an interesting conceptual aspect that the divergences in the graphs in Figs. El can only 
be renormalized consistently if the correlation of soft and ultrasoft scales, /i;/ = /i|/m = mz/^ 
is taken into account (see also Ref. 0). As an example, let us consider the contribution to 
^^NNLL,2 induced by the diagram Fig.Ek some detail. The result for the diagram reads 
{E = -{^/s - 2m + ie), = niE) 



Fig.Ek = -iCiCp 



2 OLs{l^uWs{l^s) mp 



An Att 

1 1 A. ^^Ulln^^U-- - 29 



(10) 
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-\CF{2Cl + 'iCACF + Cl) 


Cf{2CI - CaCf - C\) In 2 
-^CF{'iSCl + 27CACF + Cl) 



TABLE I: Contributions to Eq. @ from the graphs in Figs. [2 



The corresponding counterterm contribution from the graph in Fig.lHt reads 



Fig.Efc 



-iCICf 



4tt 



4tt 



(11) 



and the sum including the Ci counterterm graph from Fig.lHl gives 



Fig.|2k + Fig.Et + Fig.H 



-I C\ Cf 



Air 



1 + 12A5 



6e2 



1 /I 



e V6 



1 



p^ 



+ - ^ In^ - ^ + 2A5 In^ ln2-l -2fi5-^ + ... 



13 



1^ 



(12) 



where 6 = {C\Cf) ^- This leads to the following contribution for the Ci counterterm 

1 

~ 12 



A 
B 



C\ Cf , 



In 



7 

ln2 - - 
6 



C\ Cf ■ 



(13) 
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FIG. 4: Four-loop graphs for the calculation of the soft non-mixing contributions of the NNLL 
anomalous dimension of ci. 



(a) 




(b) 




1 T 

FIG. 5: (a) Six-field operators that contribute to the soft running of Oj^'j,-^ ^2 through the graphs 
in (b). 



The dependence on In/x^ and In fiu in B vanishes only if the correlation = l^'s/f^ is 
accounted for in Fig.|2^ as well as in the counterterm graphs needed to cancel the sub diver- 
gences. This is the case for all graphs in Figs. |21 that contribute to B. This demonstrates 
that the correlation of the soft and the ultrasoft renormalization scales, which is give n un- 
ambiguously from the mass dimension and the v counting of the operators (see Ref. 
is also needed to ensure the renormalizability and consistency of the effective theory. It also 
shows that the correlation of the soft and ultrasoft scales for all scales below m is an integral 
property of the effective theory for nonrelativistic dynamic (i.e. non-static) heavy quark 
pairs. 



B. Soft non-mixing contributions at NNLL order 

The four-loop correlator diagrams relevant for the soft non-mixing contributions are 
displayed in Figs. El where we have used again the graphical notations of Refs. [1, |^ Hs^ . 
Diagrams with insertions of two soft vertices involving the functions U'^"\ W^'^\ Y'^"'\ 
are shown in Figs.|^-e. Insertions of the 6-field operators ^^2^,2/1, 2c 

from Ref. |l3 

are shown in Fig.|3f. In Fig.E^ there are also insertions of soft 6-field sum operators (see 
Fig.E^L), which, upon closing the soft lines (Fig.EjD), contribute to the soft running of the 
operators 0]fj^^ (The ultrasoft running of the operators C^'^^ ^2 discussed in detail 
in Ref. Ts'].) The four-quark matrix elements of these operators, with the intermediate sum 
and the soft loop integration being carried out, are given in App.^ For a part of these 
matrix elements, as explained in App.El we used results obtained earlier in Ref. based 
on the asymptotic expansion of QCD diagrams. To ensure consistency of the effective the- 
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TABLE II: Contributions to Eq. from the graphs in Fig.|ll 



(1 T) 

ory under renormalization this enforces the convention for the operators (9^ ' explained 
in SeclH The diagrams needed to cancel the one-loop soft sub divergences are displayed in 
Figs.EK^g- There are no two-loop sub divergences. After the removal of the one-loop subdi- 
vergences the contribution to the renormalization factor Z^^ from each of the diagrams in 
Figs. E] can be written in the form 



47r 



A' B' 
— + — 



(14) 



The contributions to A' and B' are given in Tab.lTTl We note that the total spin operator 
for the quarks is defined as S = |(cri + (T2) using three-dimensional Pauli matrices. This is 
the common convention of non-relativistic quantum mechanics and has also been adopted 
for the vNRQCD action We also note that other schemes for the quark spin in the 

effective theory are possible. 

Like for the ultrasoft contributions, we find agreement with the prediction of the soft 
terms made in Eq. (jHJ. This demonstrates the consistency of the vNRQCD action 
under renormalization at non-trivial subleading order. In particular, the agreement shows 
that the dynamics of soft and ultrasoft degrees of freedom are both needed simultaneously 
for all scales below heavy quark mass m. 
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C. Solution of the anomalous dimension at NNLL order 



From the results shown in Tabs.|l] and |n] we find the following result for the non- mixing 
contributions of the NNLL anomalous dimension of Ci (a^ = Q;s(mz/), = ^^(mz/^)), 



NNLL, 
lci,nm ' 



2 Cf{2CI - CaCf -Cl)\n2 -\ Cf{38C^f + 21CaCf + Cl) 
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127r 



(2) 



.(2)^ 



^2 r^(2) 



TT 



2c 
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13 + y - -4 (^5 + 8S- 



4877 



+ 



101 

487r 1 ~Y 



C\ Cf + CaCI^ 



mcl \ (15) 



Let us write the solution of Eq. (jSj) for z/ < 1 as 



In 



Ci(l)J 



NNLL/ 
nm I 



(16) 



The NLL order term ^^^^{v) was determined in Refs. jl3l . Il7l |. From Eq. (|T5|) we find the 
following form for the NNLL non-mixing term, 



I- z + 2\n{w) 



+ 65 asim 



+ 66 as("^)^ 



\-2z-\z' + {A- z') Hw) 



1 _ ^2-2CA//3o 



+ 67 aJ"^) 



X _ ^2-13CA/(6/3o) 



(17) 
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where the coefficients bt read 



2Cf + Ca{S^ - 3) 



+f3oCA [ - 7ACf + Ca(42 - 13S^)] + Cl(100C^ - 9Ca)] 
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b, 
ba 



37UCaP^ 



2Cf I ^2 



-wl 



ASlCl + 350CaCf - 64C| 
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48ClC^ 



lOOCp - 9C, 



9/3o 



1 + 181n2 



QCaCf 



3 + 21n2 



-2^2 



/3o 17Cl + 30CaC7^ + 4C| 



32C^ 



19 + 181n2 
+ 3C^Ci? + 2C^ I 



C2 



144(/3o-2C^)(/3o-C^) 



3 + 4S^ 



9 - lOS^ 



+4C2 



33 + 32S^ 



Cl{hCA + 8Cf) 



39C^(6/5o-13Ca)(12/?o-13Ca) 



18/3o2 - l7lj3oCA + 385C1) 



(18) 



Note that 64 originates exclusively from the ultrasoft corrections determined in Sec. llll Al 

One can also determine the NNLL mixing contributions of the anomalous dimension at 
the hard scale u = 1 since the corresponding matching conditions of the couplings appearing 
in Eq. (jH) are available. [13, [13 The result reads 



^NNLL. 

ic\,ra \ 



1) 



487r 



Ca 16 S2 



+ ACf ( 



5-2S^ 



15t 

5 



(19) 



From Eg. p9p we can determine the as(m)^lnz/ term of the NNLL mixing contribution, 
which is the first term in the expansion in terms of a<j(m), 



er^^(-) 



7^.^(1) lnz. + 0(a^ln2z.). 



(20) 



Numerically the NNLL non-mixing contributions to ci are rather large and entirely domi- 
nated by the first term in Eq. ()15|) which originates from the ultrasoft corrections. In Tab. II I II 
the values for ^^^^{v) and ^nm'"'"(z^) are displayed for different v for the top and the bottom 
quarks. For top quarks we find that the NNLL non-mixing contributions are of the same 
size as the NLL terms for the relevant region v v 0.2. Here, the new NNLL order 
corrections shift Ci by about -|-5%. However, it is not yet possible to draw phenomenological 
conclusions for the normalization of the top threshold cross section in e^e~ collisions from 
this result, since the NNLL mixing contributions for v < 1 are still unknown. Because 
these corrections involve two ultrasoft loop integrations arising in the NLL running of the 
coefficients in Eq. (jU, they could potentially be large as well. It is therefore an important 
task to determine the mixing corrections for all scales below m. 

For bottom quarks the NNLL non-mixing contributions are more than five times larger 
than the NLL terms for the relevant region z/ ~ w ^ 0.3-0.4. This is not unexpected 
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because for bottomonium systems the binding energy ~ mv"^ is already of order Aqcd- Thus 
our result seems to affirm that for bb states non-perturbative effects have a rather strong 
influence, and that the vNRQCD description ceases to work even for the ground state. jl3| 
However, also for the bottom quark case the knowledge of the NNLL mixing contributions 
might be useful to gain further insight into this issue. 





m = 175 GeV 


m = 4.8 GeV 












1.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.9 


0.0033 


0.0019 


0.0145 


0.0181 


0.8 


0.0069 


0.0041 


0.0308 


0.0425 


0.7 


0.0110 


0.0069 


0.0495 


0.0770 


0.6 


0.0157 


0.0104 


0.0712 


0.1304 


0.5 


0.0211 


0.0150 


0.0968 


0.2254 


0.4 


0.0274 


0.0216 


0.1335 


0.4537 


0.3 


0.0349 


0.0318 






0.2 


0.0435 


0.0512 







TABLE III: Numerical values for ^'^'^'^(z^) and ^nm^^(z^). The values for m are pole masses. The 
numbers are obtained by evaluation of the analytic results using four-loop running for as and 
taking ai^^ (175 GeV) = 0.107 and ai''^ (4.8 GeV) = 0.216 as input. 



D. Production Cross Section at Threshold and Comparison 

Due to the convention we use for the order a^f ° potentials, our formulas for the NNLL 
order vector current correlator Ai, which is used to express the quark pair production 
cross section in e"'"e~ annihilation, slightly differs from Ref. This also affects the two- 
loop matching condition for ci. The difference arises because the corresponding matrix 
elements are UV-divergent and the d-dependent contributions that arise from summing the 
intermediate indices in the operators Ol''^ in dimensional regularization (see App.^ lead 
to modifications of the UV-finite terms. Altogether, we now need four different types of 
corrections to the current correlator, SG^j^q^, ^G^pg, 6G''^ and 6G'''^ to account for the 
corrections originating from the order a^v^ potentials. The results for 6G^^ and were 
given in Ref. [13[. The corrections SGqj^^q^ and 5G'cf2 arise from the terms proportional to 
GaCf and Cp, respectively, in the color singlet combination of the operators Ol and 
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and have the form 



m 



+ 



57r a 

2 

m 

Svr a 

2 

m 



+ 



57r a 

2 

m 

Stt a 



IV — a 



IV — a 



In 



—2 1; 



/I 3 



16 Ve 



In 



5 ( % a 
- + ln2 + 7£; + ^ 1 
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-« V 



1 + In2 + 7i5 + ^ 1 



2^ 



a' (\ 2 



16 Ve 



+ T77 12 



The vector current correlator A\ at NNLL order then reads 

Ax{v,m,v) = 6iVe G\c;,v,m,v)^(v^^\v)^2V^;\v)^bG\a,v,m,v) 
+ V^") (z/) 5G"(a, V, m, z/) + 5G'^^^(a, v, m, z/) 



— CaCf ag{miy) 6G(^js^Qp{a, v, m, z/) + Q;^(mz/) ^Gcpgl*^) ^? ""^j 
+ a2(„ii.)vf^)(z/) 6G''\a, v, m, v) + al{mv)V^^l{v) bG^''{a, v, m, v) 
a = — —Vl'Av) 



4tt 



47r 



2Cf 



ln2 31 2 

~ ~ 24 ~ ^ 



(21) 



(22) 



where 56"^, (56"-, ^G'^'" and (5G"=^ vi'\ were given in Refs. and 0, 

respectively. The two-loop matching condition for c\ now reads 



, (23) 



where the constant k was determined in Refs. |1S] . 

Recently, the order a^lna^ corrections to the heavy quarkonium partial width into a 
lept on p air were computed using an asymptotic expansion of QCD diagrams close to thresh- 
old [19^. In this work the summation of higher order logarithms was not attempted. The 
results were partly based on three-loop quark-antiquark-to-vacuum matrix elements in the 
ofT-shell limit mE 7^ 0, = |2^. With the result for the total cross section of quark 
pair production at threshold in e^e~ annihilation in Refs. 0,0,^^! and including the new 
formulas given above, the order o?^ Ina^ corrections can also be derived by expanding out the 
summations contained in the Wilson coefficients for energies on the bound state poles. The 
contributions proportional to = S'(S' + 1) in our result do not agree with Ref. |l9| because 
there the terms S'(S' + 1) do not refer to the three-dimensional quark spin operators ji^ and 
are defined in a different scheme for the quark spin. For the In corrections to the heavy 
quarkonium partial width into a lepton pair we find a discrepancy, which would correspond 
to the term \C\Cf - \CaGI{2 - 5 In 2) being added to the constant Ci of Ref. H^- This 
term has the same sign and approximately the same size as the constant C\ itself and affects 
the error estimates made in Ref. [l9l |. 



14 



IV. CONCLUSION 



In this work we have determined, within the framework of vNRQCD, the non-mixing 
contributions to the NNLL anomalous dimension of the leading spin-triplet current Ji,p, 
which describes the production of non-relativistic quark-antiquark pair in e~^e~ annihilation. 
Our result for the terms of the three-loop renormalization constant of the current is 
consistent with the constraints of renormalization group invariance imposed on the already 
known NLL UV-divergences. This demonstrates that both soft and ultrasoft degrees of 
freedom need to be present in the effective theory for all scales below the heavy quark 
mass m as predicted by vNRQCD. Moreover, a sequence of different effective theories with 
an intermediate uncorrelated matching scale is not consistent under renormalization. Our 
results also show that the NNLL order renormalization constant of the current is independent 
of the renormalization scales for the soft and ultrasoft degrees of freedom, only if their 
correlation is accounted for, as given unambiguously from the mass dimension and the v 
counting of the operators in vNRQCD. We have derived an updated expression for the 
NNLL total cross section of heavy quark pair production at threshold in e"'"e~ annihilation. 
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APPENDIX A: CONVENTION FOR THE POTENTIAL OPERATORS O 



The order a^v^ potential sum operators O^. ' used in this work are obtained by matching 
to the one-loop QQ scattering amplitude in full QCD in the threshold expansion and carrying 
out only the dq^ integration of the loop momentum q = {q^-, q) for the contributions of order 
a^f". The explicit expression for the operators is [qs = gsifriv)) 



O 



(1) 



4 ~4e 

Am 



5'o - 92 



[^i'^px-p'X-p] 



p,p ,q 



O 



(T) 



^ r l,„ „ , „ 1 

4m 



- I^^'^ - ^0 + Ca9i + -AC, + Ca) 92 [^^,T^^pxLp.T-^X-p] 



p,p ,q 



(Al) 



where the functions Qi have the form 



^0 



9i 



92 



[q-pj 



iq-p 



l\2 



IP-P 



i\2 



(q - p)' + (q - p Q 
(q-p)Mq-pO 



l\2 



(q-p)' + (q-pT 



IP-P 



l\2 



(q - p)^ (q - p')^ (p - p') 



2 2 

— 



+ 



q2 _ p'2 



(q_p)4(q_p,)2 (q-p)2(q-p') 



(A2) 
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They give the contribution A^p = ^C^^^ + V^pof^ to the vNRQCD Lagrangian. The 
Wilson coefficients read 



vfV) 



1. 



(A3) 



Note that the tilde over the renormalization scales iiu and iis refers to the MS definition 



/i[7 5, where p = — ln(47r). When O), ' are inserted in four-quark matrix 
elements or used in the Schrodinger equation with the sum over the intermediate index q 
being carried out in dimensional regularization, they take at leading order the form 



4 ~4e 



4m(k2)(2-«/2) 
4m(k2)(2-W2) 



/(1,1)-/(2,1) 



I® I. 



-Ca - 2Cf /(1, 1) -[-Ca- 2Cf /(2, 1 



(A4) 



where 



/(a, 6) 



r(a + &-f)r(f -a)r(f -6) 

r(a)r(6)r(n-a-6)(47r)t 



(A5) 



and k = p — p'. Here n = (i — 1 = 3 — 2e. The matrix element of the color singlet combination 
O 



is) 



Ol}^ - CfO^ reads 



m 



(^k2)(2-n/2) 



— CaCf -\ 



(A6) 



which agrees with Ref. |15| . 



APPENDIX B: 4-QUARK MATRIX ELEMENTS OF 6-FIELD OPERATORS 

(1 T) (1) (T) 

The soft running of the operators Of, , O)^^ and O)^^ originates from four-quark 
matrix elements of 6-field sum operators (see Fig.|3^) in close analogy to the operators 
0^^2A2c-i "which contribute to the soft running of the spin-independent order agV (l/m^, 
(p2 + p'2)/(2m^k2), S"^ I'm?) potentials. j3l The divergences in these 4-quark matrix ele- 
ments have already been deduced earlier, but for the computation of the NNLL anomalous 
dimension of Ci their full form is required. To be definite we call the 6-field operators 
Ol' , Of,_^ and O^J. The operators ' and , Oj,J are responsible for the soft 

(1 T) (1) (T) 

running of the operators Ol' and Of,_^ , ' ^respectively. They give the contribution 
ACp = Vi^^Ol}^ + Vr^T + ^ki^kl + W2^T to the vNRQCD Lagrangian. The op- 

~ (1) ~ (T) 

erators 0\. and 0\ arise from a matching computation at two loops similar to the one 
described in App.lXlfor In analogy to Eq. ()A3|) one finds that VP{h') = Vj^\h') = 1. 

The leading order 4-quark matrix element for the unrenormalized color-singlet combination 
of C^''' > with the sum over intermediate indices and the soft loop being carried out (Fig.|3|D) 
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can be derived in an expansion in e from results given Ref. M and reads 



7ra,(mz/)3/2|^/if M 
m(k2)7/2-n P \ e 



Ca[ — + ln2 I +Cf\- - -ln2 
^'48 / V3 2 



15 4 

16 ~ 3 



-Cl ( - T= ln2 ) - CaCf ( 3 - ^ ln2 



(Bl) 



which is sufficient for the determination of the renormahzation constant for Ci. The divergent 
/3o term is responsible for the running of as contained in the definition of O^^^ (Eq. ()A6|) ). 
while the other divergence is related to the evolution of V^/ and V^2 • 

The renormahzation of and 0)^J is in complete analogy to the one of the operators 

O^P^ and Of^. Thus one finds V^lV) = "^kii^) and V^a V) = V^P(z/) at LL order, and the 
coefficients vanish at the hard scale. The leading order four quark matrix elements have the 
form 



26-n, ■ al{mu)Ti^lL%' (1 - nf (12 - 6n + n^) 



m(k2)'^/2^" n{n — 2) cos 



2 ' 



(4r2-l)- 4 Triz 



4_„ . al{mu)Ti^lL%^ (1 - n)2 (17 - 9n + 2^^) 



m(k2)7/2- 



n(n - 2) cos(^) 



CA(4n-l)-4Tn; 



Al)/f^,l/2)T- 



(B2) 



The divergences are responsible for the running of the coefficient Vc contained in the 
and Og^ 



definition of CiV and OZ' . 



APPENDIX C: COLLECTION OF WILSON COEFFICIENTS 

The gauge invariant HQET coefficients that appear in this work are 3 



-Ca/Po 



13 13C 



A 



-13Ca/(6/3o) 



(CI) 
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For the color singlet channel the coefficients of the order agV potentials relevant for our 
result are [§, [l^ 



1 1 ^ ^ 

1 :r-^ m[Wj 



V^'V) = 7TCFas{m) (z-l) 



33 32Cf 9Ca IOOCf 



13 13Ca 13/?o 13/?o 
87rCj7(3/3o - 1 ICa) {5Ca + 8Cf 



13Ca(6/3o-13Ca) 



asim 



7iCf(J3o^a) 
(Po-^Ca) 

-2ttCf 



ajm z 



^^1-2Ca//3o _ ]_j 



^^l-(13CA)/m) _ ij 

1QttCf{Ca-2Cf) 
Wo 



as{m) zln(w) 



as[m) 



1 



{2Ca-Po) 

whereas for the order a^f ° potentials they read 



CA+-{2f5,-lCA) 



(C2) 



8CaC 



- \n{w) 



^^(^)^ \ 2CACFiCA + Cd) , - , . 

n2\^) = H^)- (C3) 



The coefficients associated to the 6-field operators 



(2) 

2<fi,2A,2c 



3/?o 
read |l 



3Po 

-^Ca 
3/9o 



<^26 l'^) = ^ Hw) 



3Po 



ln(w) 



For the results we used the definitions 

as{m) 



w 



(C4) 



(C5) 



For SU(A'"c)-QCD the color coefficients that appear are 

— I I 1 2 

Ca = Nc, Cf = — , T=-, Cd = 8CF — SCA, Ci = -CfCa~Cp. (C6) 
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